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Abstract. We re-examine the importance of slow- roll corrections during the evolution 
of cosmological perturbations in models of multi-field inflation. We find that in many 
instances the presence of light degrees of freedom leads to situations in which next 
to leading order slow-roll corrections become significant. Examples where we expect 
such corrections to be crucial include models in which modes exit the Hubble radius 
while the inflationary trajectory undergoes an abrupt turn in field space, or during a 
j> phase transition. We illustrate this with several examples - hybrid inflation, double 

quadratic inflation and double quartic inflation. Utilizing both analytic estimates and 
full numerical results, we find that corrections can be as large as 20%. Our results 
have implications for many existing models in the literature, as these corrections must 
be included to obtain accurate observational predictions - particularly given the level 
of accuracy expected from CMB experiments such as Planck. 

Keywords: inflation, cosmology of the very early universe, cosmological perturbation 
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1 Introduction 



The inflationary paradigm [1,2] remains a convincing causal mechanism for providing 
the needed initial conditions of the early universe - despite scrutiny from a wealth of 
precision cosmological observations [3, 4] (see [5, 6] for reviews). Inflation accomplishes 
this by providing a period of quasi-de Sitter expansion that leads to a classical spec- 
trum of nearly scale invariant density fluctuations [7]. The evolving density contrast 
then leads to the large scale structure and Cosmic Microwave Background (CMB) 
anisotropies observed today. In the simplest models of inflation the acceleration is 
assumed to be the result of a single scalar degree of freedom - the inflaton. In such 
models it was shown long ago that the perturbation in the spatial curvature is a con- 
served quantity on super-Hubble scales [8] , which makes single field models of inflation 
rather easy to analyze. 

Embedding the concept of inflation in a more fundamental theory generically leads 
to the existence of additional light degrees of freedom. These fields typically influence 
the dynamics and may even help drive inflation, leading to multi-field models (early 
examples are [9, 10]). In such situations the curvature perturbation does not necessarily 
remain constant on super-Hubble scales, resulting in new theoretical challenges as well 
as richer possibilities for observations [11-19]. One example arises when the inflationary 
trajectory deviates from a geodesic in the field space [20, 21]. As a result, the curvature 
perturbations are not constant after crossing the Hubble radius, and the naive single- 
field picture fails. Thus, the usual single-field relations between observable quantities 
(such as the amplitude of the perturbations and the spectral index inferred from the 
two-point correlation function) and the parameters of the inflationary potential are 
lost. Another example is when a time dependent coupling between fields approaches 
the regime of strong coupling, which in turn can lead to interesting observational 
signatures [20-25]. 

Among the theoretical obstacles facing multi-field inflation models is finding sim- 
ple (and universally applicable) tools for calculating the power spectrum of curvature 
perturbations, particularly in the situation where couplings between fields are signif- 
icant. A common and useful tool for analyzing multi-field models is the so-called 
SN formalism [26-28], which allows to calculate the power spectrum of the curva- 
ture perturbations to leading order in the gradient expansion and to all orders in the 
slow-roll parameters. Indeed, this approach has led to interesting and quite restrictive 
bounds on multi-field inflation models from primordial non-gaussianity [29-32]. How- 
ever, using this method in practice requires (typically reasonable) assumptions about 
the amplitude of the fluctuations at Hubble radius crossing. In particular, one assumes 
that the amplitude is equal to that of perturbations of a light scalar field in pure de 
Sitter space, where leading order in slow-roll is more than adequate to capture the 
dynamics - an intuitive and very accurate approximation in many cases. However, in 
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special cases of hybrid inflation this approximation may fail, specifically when mildly 
tachyonic isocurvature modes source the curvature perturbations [33-35]. These spe- 
cial situations provide examples where the leading slow-roll approximation can fail to 
provide accurate initial conditions for the SN formalism. 

With these motivations in mind we revisit multi-field models of inflation. We find 
the next to leading order corrections to the power spectra of curvature and isocurvature 
perturbations near Hubble radius crossing and on super-Hubble scales utilizing the 
transfer matrix formalism. We give simple examples of models where these corrections 
are essential for obtaining accurate estimates of the power spectra. This generalizes the 
work of Refs. [36-40] to the multi-field case, and extends the results of Refs. [16, 41] 
to include higher order corrections and an explicit expression for the transfer matrix 
of perturbations. 

Outline. — The paper is organized as follows. In §2 we review canonically normalised 
two-field inflation, discuss the usual slow-roll approximations and motivate considering 
higher order slow-roll effects. In §3 we study the dynamics of the perturbations on a 
Friedmann-Robertson- Walker background, around Hubble radius crossing. We derive 
the expressions for the power spectra of curvature and isocurvature perturbations up 
to next-order in the relevant slow-roll parameters. These analytical results are valid 
up to a few e- folds after horizon crossing (when the modes have become classical), and 
capture the relevant time and scale dependence of the power spectra. The evolution of 
the perturbations on super-Hubble scales is the subject of §4, where we calculate the 
transfer matrix for the perturbations up to next-next-order in slow-roll parameters. In 
§5 we present two simple but instructive examples of inflationary trajectories in hybrid 
inflation and double quadratic inflation models, for which the importance of the next- 
order corrections is manifest. We conclude in §6. The appendix collects additional 
formulae necessary to compute the power spectra. 

2 Perturbations in multi-field inflation 

We will restrict our attention to the case of two scalar fields minimally coupled to 
gravity, although our analysis is easily generalizable to more fields. At the level of 
renormalizable interactions the action is then 



where m p = 1 /\/8itG ~ 2.43 x 10 18 GeV is the reduced Planck mass. The background 
equations of motion for the fields = {<p, x} are given by 
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where the dot denotes differentiation with respect to cosmic time and = dV/d<fi. 
The cosmological evolution is described by the Friedmann equations 

W*ml = <?+x i + V(<p,x), (2-3a) 
2Hm 2 p = -[ip 2 + x 2 ]i (2.3b) 

where H = a/ a is the Hubble parameter and in what follows we will work in units 
where m p = 1 for simplicity. We now expand the fields in fluctuations around their 
homogeneous background values 

y>(t,x) = <p (t) +6<p(t,x), (2.4a) 
x(t,x) = xo(t) + 5x(*,x). (2.4b) 

A particularly convenient basis for interpreting the behavior of the fluctuations [12] is 
found by performing the instantaneous rotation 

5a = cos 9 5ip + sin 9 5x , (2.5a) 
5s = — sin 9 5ip + cos 9 5x ■ (2.5b) 

The rotation angle 9 is given by cos 9 = tp/a and sin (9 = x/<7, where a = \/(p 2 + x 2 - 
The background equations of motion (2.2) then become 

a + 3Ha + V a = , (2.6a) 

9 = , (2.6b) 
a 

where V a = cos 9 + sin 9 V x and V s = — sin 9 V jip + cos 9 V x } From the background 
equations we see that in this basis 5a corresponds to the scalar field fluctuations along 
the trajectory of the inflaton (adiabatic perturbations), whereas 5s is the fluctuation 
orthogonal to the trajectory (isocurvature perturbations). 

Whereas the entropy perturbation 5s is automatically gauge invariant, we need 
to introduce a gauge invariant definition for the adiabatic perturbations. This may be 
accomplished by introducing the Mukhanov-Sasaki [42, 43] variable Q a = 5a — 
which corresponds to the instantaneous curvature perturbation on surfaces of constant 
a, where $ is the Bardeen potential describing scalar perturbations of the metric 2 . 

Given the gauge invariant perturbations Q a and 5s, it is convenient to introduce 
the rescaled variables u a = clQu and u s = a 5s, in terms of which the equations of 
motion for modes with comoving momentum k take the form [17] 



a 



a" 



+ a 2 C a 



u a + 



2a' -± + a 2 C a 
a 



u s = (2.7a) 



1 More generally, with I,j£ {a, s} and i, j £ {ip, x}, we have Vjj,,,k = E}Ej . . . Ej ( V^j...k, where 
= E} = cos (9, E% = -Ef = sin 6. 

2 On super-Hubble scales, Q a is related to the comoving curvature perturbation via 1Z = (H/&)Q a . 
Similarly, we can define a comoving isocurvature perturbation through S = (H/a)Ss. 
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and 



- 2a— u„ + 



a 



a" 



k + a z C s 



u s + 



2d + a 2 C s 
a 



u a = , (2.7b) 



where primed quantities are differentiated with respect to conformal time r = J dt/a(t) 
and the coefficients Cu are listed in Eqs.(A.l). We note that, since we are restricting 
our attention to classical scalar fields with renormalizable interactions, both the scalar 
fluctuations u a and u s propagate at the speed of sound c 2 = 1. 

2.1 Slow-roll approximation: the need for next-order contributions 

We can simplify the equations of motion for the perturbations (2.7a) and (2.7b) by 
applying slow-roll conditions. In many realistic models of inflation, the approximate 
constancy of the Hubble parameter requires e = —H/H 2 <C 1 , while the approximate 
scale invariance of the power spectrum of perturbations at Hubble radius crossing 
(proportional to H 2 /e) requires the smallness of 

^ = 2(25-^), (2.8) 

where rju = Vjj/3H 2 . This yields e, \r) aa \ <C 1. These requirements, however, do not 
place severe constraints on the remaining slow-roll parameters, since the time evolution 
of t)h is given by [16]: 

^■q aa = 2(67]^ - r]l s ) - Caa , ( 2 -9a) 

jjV™ = V*s{2e + r) aa - rj ss ) - g as , (2.9b) 

^Vss = 2(e Vss + nl,) - Css , (2.9c) 

where $ JK = V^Vjjk /V 2 . The quantities 3H 2 r]u contribute to the mass matrix of 
the perturbations. In particular, rj as parametrizes the coupling between curvature and 
isocurvature perturbations and it is therefore related to the bending of the inflationary 
trajectory in field space via 9/H = —rj as /{\ — e/3). 

As also observed in [25], there are three patterns for the slow- roll parameters e 
and rjij consistent with the requirements mentioned above. One possibility is that 
they all have magnitudes much smaller than one. This is the pattern most often 
considered in the discussion of two-field inflationary models. Another possibility is 
that e, \T] aa \ ^ 1, rj as ~ and r] ss is positive but otherwise arbitrarily large. This 
is the decoupling regime, which is effectively equivalent to a single-field inflationary 
model. A pattern interpolating between these two possibilities is also conceivable with 

(£, \V<ra\, \r]*s\) < \Vss\ < 1. 

If all slow-roll parameters are much smaller than one, the evolution of the modes 
is simplified: the perturbations freeze-in near Hubble radius crossing, and the effect of 
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a small and almost constant coupling between curvature and isocurvature modes can 
be easily included [15]. 

A quantitative description becomes more intricate if there is a significant super- 
Hubble evolution of the isocurvature modes, as in the recently considered models of 
hybrid inflation [33-35], allowing for several tens of e-folds of inflation after a mild 
tachyonic instability develops, associated with a moderate (and negative) i] ss . One 
can, however, get some insight into the inflationary dynamics by solving Eqs. (2.9a)- 
(2.9c) in the formal limit of £ — > and with negligible With initial conditions 

r\ aa 0, rips ~ and f] ss = rjis < at early times, the solution is: 



I Iss 

V&cr,ss ~~ ~ 



lTtanh[^?(iV-JVo)] 



„(o) 

Vas = ± , (2.10) 

2cosh [^(N-No 1 



where N is the number of e-folds and N is a constant of integration. Hence, although 
the hierarchy (e, \r) ar7 \, \T) as \) <C \r] ss \ < 1 is not stable if T] ss < 0, it can persist for 
a large number of e-folds and the transition to a stable one with 7] aa = r]is* takes 
approximately l/|r/i^| e-folds. It is therefore interesting to study the inflationary 
dynamics with moderate and negative t] ss . 

Another instance in which the dynamics can be quite involved is the case when the 
slow-roll parameters are small only initially, with some becoming of order one at a later 
time, without violating the slow-roll condition e < 1. This is the case, for example, in 
the model of double quadratic inflation [44] or double quartic inflation. With rj aa = cue, 
as it is the case for a single monomial dominating the inflationary potential, we can 
initially neglect rj as , r] ss , ££ ffa and £% ss in Eqs. (2.9a)-(2.9c), obtaining 

l-2e(°)(2-a)(JV-iVo)' V " Vas \eWJ ' 1 ' 

which means that a sudden growth in e and r] aa drives a subsequent growth of rj as and 
r] ss until these slow-roll parameters backreact on e and rj a(T and the approximation used 
here breaks down. 

We emphasize that in all of these models, in order to describe accurately the 
evolution of the modes outside the horizon, one must keep higher order terms in the 
slow-roll parameters. This will be the focus of §5. Motivated by the discussion above, 
our strategy for probing the dynamics of the curvature and isocurvature modes will 
be the following. First, in §3 we will track the evolution of the perturbations up to a 
few e-folds around Hubble radius crossing. At this stage all slow-roll parameters are 
restricted to have magnitudes smaller than one, with the notable exception of r] ss . As a 
result, the solutions of the equations of motion (2.7a)-(2.7b) will have to be expanded 
to next-order in rj ss , but only to leading-order in the remaining slow-roll parameters. 
In §4 we will then follow the evolution of the modes outside the Hubble radius. Since 
in this region the term proportional to k 2 in Eqs.(2.7a)-(2.7b) becomes negligible, the 
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solutions simplify significantly and enter the non-oscillatory regime. However, as we 
will see, to obtain a faithful estimate of the amplitude of the perturbations one generally 
needs to expand the equations of motion to next-order in all slow-roll parameters. 



3 Dynamics near Hubble radius crossing 

We start by noting that the expansions to next-order in slow-roll of the Cjj coefficients 
in Eqs.(A.l) and of the background scale factor in Eqs.(2.7a)-(2.7b) do not contain any 
terms proportional to rf as . As a result, these quantities can be expanded to leading- 
order in slow-roll. The first order derivative terms in Eqs.(2.7a) and (2.7b) can be 
removed by changing variables to 



w\\ ( cos $ u a + sin i? u s 
w 2 ) I — sin d u a + cos u s 



(3.1) 



where $ = —V s /&. As we justify in the Appendix, the angle •& can be chosen so that 
the equations of motion for w\ and w 2 decouple when a given reference scale, k*, crosses 
the Hubble radius 3 . The equations of motion around the time of horizon crossing for 
k* then simplify to 



// 

2 



W 



+ 



r 2 V 2 + A 2 * 



w 2 



, (3.2) 



where Ai* and A2* are functions of the slow-roll parameters given in the Appendix. 
3.1 Slow-roll expansion of the mode functions 

The equations of motion (3.1) for the variables wa encapsulate the dynamics of the 
curvature and isocurvature perturbations around Hubble radius crossing. The positive 
frequency solutions to the equations of motion (3.2) are represented, up to an irrelevant 
phase, by 

WACx y^^H^e A (k), (3.3) 

where e^(k) are Gaussian random variables 4 such that (e J 4(k)e^(k / )) = Sab S^(k— k'). 
Here A,B e {1,2} and Hj$ are Hankel functions of the first kind and of order 



,H = X I^ + 3\ A ~1 + \ A -l\* A . (3.4) 



3 Here and throughout the paper quantities with the subscript * are to be evaluated when the 
reference mode crosses the Hubble radius, i.e. when /c* = a^iJ*. 

4 Since our goal is to calculate two-point correlation functions at tree level, the use of classical 
random fields here is equivalent to the proper procedure of expressing the Fourier modes in terms of 
creation and annihilation operators. Such a representation is legitimate on super-Hubble scales, so 
this choice will allow us to make an easy comparison with the results of §4. 
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We emphasize that the presence of \\ in the expansion (3.4) introduces terms propor- 
tional to if ss in the expansion of the Hankel function in the slow-roll parameters. We 
also note that deep inside the horizon, when k ^> aH, the modes are rapidly oscil- 
lating and the solution (3.3) reduces to the standard Bunch-Davies vacuum [45]. The 
two-point correlations of the w A variables satisfy 

(w\w B ) = - A {-T)\H^{-kr)\ 2 5 AB , (3.5) 

in terms of which we can write the correlation functions for the u a and u s perturbations 
[36]. 

Our goal is to obtain approximate expressions for the correlation functions in 
Eq.(3.5) around Hubble radius crossing. This exercise is particularly simple in single- 
field inflation models where the comoving curvature perturbation is conserved on super- 
horizon scales. There, the limit \kr\ — > gives both the correlation function value a 
few e-folds after horizon crossing and its asymptotic limit (with only the growing 
mode left) . Indeed, if we take the decaying mode into account 3 to 5 e-folds after 
horizon crossing, the asymptotic limit of the Hankel function (with order approximately 
3/2) will be corrected by O(10~ 6 — 10~ 4 ) terms. In contrast, if inflation is driven by 
multiple fields some of the slow- roll parameters might be non- negligible, which allows 
a significant deviation of the order of the Hankel function from 3/2, which induces a 
much stronger time dependence than the presence of a decaying mode. In order to 
estimate this effect, we expand the Hankel function up to second order in the slow-roll 
parameters, keeping only the constant terms and the dominant contributions when 
\kr\ -»■ 0. We find that 

H^U-kr) | 2 .- 1 1 - 2Ayu [ - 2 + 7E + ln(-2*r)] + X 2 A J(-kr) } , (3.6) 

where 7e — 0.577 is the Euler-Mascheroni constant and the time dependent function 
f(—kr) satisfies 

6/(-fcr) = 16 + 37r 2 + 4( 7E + ln2)(-ll + 3(7E + ln2)) + 
+ (4 + 6(7e - In 2)) ln(-frr) + 12 ln 2 (-fcr) . 

The discarded quadratic contributions in slow-roll parameters above are of the same 
order as the discarded decaying mode. Therefore, this asymptotic expansion is valid 
up to 1 part in 10 4 and, as we will see in §4, it encodes the relevant dynamics needed 
for an accurate description of the evolution of the modes outside the horizon. 

Although the asymptotic expansion of the Hankel function in Eq.(3.6) is often a 
fairly accurate approximation of the predictions obtained by solving the full equations 
of motion (2.7a)-(2.7b) for the curvature and the isocurvature perturbations, it has 
some limitations. One of these comes from the expansion of the scale factor 

1 



2 / (3-7) 



a(r) 



1 + e± + e± ln(-fc*r) 



(3.8) 
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which is only valid for | ln(— k±r)\ < Secondly, the diagonal form of the mass 

matrix in Eq.(3.2) is strictly true at Hubble radius crossing. Since this matrix is rotated 
by an angle ~ rj as ln(— k+r), our result is only limited to times for which this angle can 
be considered small, unless Ai* ~ A2*. Finally, as the slow-roll parameters evolve slowly 
in time, we can only treat them as constants when ln(— k+r) is much smaller than any 
of the right hand sides of Eqs.(2.9a)-(2.9c). 

3.2 Power spectra 

The two-point correlation function of a given perturbation Q is related to the (dimen- 
sionless) power spectrum via 

2tt 2 

(QaQl) Q B (k')> = (27r) 3 5( 3 )(k + k')—V AB (\k\) , (3.9) 

where Qa,b = 11, S (repeated indices are customarily shortened to a single one). Eval- 
uating the power spectra around the time of horizon crossing, we find 5 

H 2 ( 1 

V n = < 1 - 2e * + 2 [2 - 7E - In 2 + \n(k*/k)] [3e* - W ] \ , (3. 10a) 

H 2 f 1 

Vns = - 2 + 7e + In 2 + In(-kr) | V(TSi! , (3.10b) 

H 2 ( 

V s = mt- i 1 - 2e * + 2 H-hr) [2e* - W + r) ssi ] 

+ 2(8, - r] ssir )(2 - 7E - In 2 + Hh/k)) + vLfi-kr)) , (3.10c) 



where V-r, denotes the power spectrum of the comoving curvature perturbation, Vs 
that of the isocurvature perturbation, and Vns is the cross-correlation between the 
two perturbations. 

Note that there are two types of logarithmic contributions in the power spectra: 
those which depend explicitly on time, of the form ln(— kr), and those which are scale 
dependent, ~ ln^/fc). 6 We can eliminate the latter type of logarithmic terms, by 
choosing an appropriate reference scale, in this case k+ = k. One still has to deal, 
however, with the time-dependent logs, which cannot be removed. Although it follows 
from Eq.(3.8) that at Hubble radius crossing — kr = 1 + 0(e+) and the time-dependent 
logarithms vanish (up to corrections of higher order than we are considering here), we 
will see in §5.1 that their inclusion beyond leading-order may be essential for accurately 
tracing the evolution of the power spectra for a few e-folds after Hubble radius crossing. 

Finally, the absence of time-dependent terms in Eq.(3.10a) is just a reflection 
of the fact that, without isocurvature modes, the curvature perturbations are frozen 



5 This clarifies previous results found in Refs. [16, 17]. 

6 The appearance of both families of logarithmic contributions was mentioned in Ref. [40], whereas 
Refs. [37, 46, 47] obtained time dependent logarithms. 
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after Hubble radius crossing. We emphasize that the isocurvature modes will source 
the curvature perturbations, but this is a next-order effect in the slow-roll parameters 
not involving rj ss and can therefore be neglected around Hubble radius crossing, as we 
have argued previously 7 . However, when describing the entire inflationary period this 
sourcing cannot be neglected, as we will see explicitly in the next section. 

4 Dynamics after Hubble radius crossing 

The results derived in §3 are valid up to a few e-folds after the relevant scales have 
exited the Hubble radius and the comoving curvature perturbation 1Z has become 
classical. To follow the subsequent evolution of the modes, however, one needs to 
resort to other techniques such as the 5N formalism or the transfer matrix method. 

To track the evolution of the perturbations after Hubble radius crossing it is 
convenient to write the equations in terms of the Mukhanov- Sasaki variables Q a and 

5 s [17]: 

(t)^|S)(tH^fe £)] ft)-- («> 

The coefficients Cu are given by Eqs.(A.l), and can be expanded consistently to 
next-order (see Appendix). On super-Hubble scales we can neglect the k 2 /a 2 term 
in Eq.(4.1). Furthermore, noting that in this regime the variables Q a and 5s change 
slowly compared to the scale factor, we can also neglect the double time derivatives. 
Thus, using Eqs.(A.2a)-(A.2d) we obtain, to leading-order in the slow-roll parameters: 



Qa I H = - (rj a<7 - 2e)Q a - 2 Vas 5 s , (4.2a) 
5s /H = - Vss 5s . (4.2b) 

Differentiating Eqs.(4.2a)-(4.2b) with respect to time and using Eqs.(2.9a)-(2.9c), we 
can obtain next-order slow-roll expressions for Q a and 5s. Plugging them back into 
Eq.(4.1) and again using Eqs.(A.2a)-(A.2d), we find a next-next- order slow-roll ap- 
proximation for the full equations of motion: 

AQ a + B5s, (4.3a) 
D5s, (4.3b) 

7 We expect Vm to evolve as the inflationary trajectory bends in field space. In particular, since 
9 ~ —Hrj as , its logarithmic time dependence will be proportional to rfe s . 



d Q a 

dN 
d5s 

dN 



-9- 



where we have defined 

A = -q aa + 2e- l -rf aa + ^£ Vaa - ^e 2 - ^ s - l -£ aa + 6&A , (4.4a) 

2 2 2 

B = -2 Vas + 2er] as - -r) as r] aa - -r) as r) ss - -& aB + 5 {3) B , (4.4b) 

D = -Vss - - \vs S + \eVss ~ \ilss + 8®D , (4.4c) 

where 5^ A, 5^ 3 'B and 5^D denote next-to-next- order corrections given by: 



95^ A =28e 3 - 49e V + 23^ - 2 V 3 aa + I5erj 3 as - A Vaa rjl s - 2^ M + 

+ 4eg aff - ^£ au - 2 VrTS Ca S + 3^ ff , (4.5a) 
95 {3) B = - 58e 2 r] as + iOeT] aa rj ss - Arf a(T r) cs - Arf as + 2Aer] as r] ss - ir] acT rj as rj ss - 

+ 2SL S , (4.5b) 
95 {3) D =Wer)l s - 2rj aa r]l s - 5e 2 r] ss + 2er) a(J r) ss - 4r/ CTS r/ ss + 5er] 2 s - 2rf ss - 

- QvULs ~ ^sCss + S L , (4.5c) 



where S^ 7J = ^ 2 U /H. The next-order result has been first obtained in Ref. [41]. Our 
inclusion of the next-next-order result allows us to estimate the accuracy of the next- 
order calculation and to discuss general conditions for the applicability of the technique 
described here. 

Eqs.(4.3a)-(4.3b) can be integrated to give 



Qa{N) 



N 



Q a * + 5s, / dN"B(N")e^ ^(d(n')-A(n')) 



e JlwA(N )f (46a) 



5s(N) = 5s*e^ A * D ^\ (4.6b) 

where Q aic and 5s* are the initial conditions at N±. Taking N* to be the e-fold at which 
the mode of interest crossed the Hubble radius and applying Eqs.(3.10a)-(3.10c), we 
find that the instantaneous power spectra for 1Z and S can be written as 



■N „ \ 2 pN „ 

!+(/ dN"B(N")e& dN ' GiN,) ) +cj dN"B(N")e& dN ' G{N,) 



(4.7a) 

V s = Vs,e^l AN ' G ^'\ (4.7b) 

where G = D — A, V* = H 2 /87t 2 e+ and Vs* is given by Eq.(3.10c) setting k+ = k and 
— kr = 1. The relative correlation between curvature and isocurvature perturbations 
at Hubble radius crossing is denoted by c+ = —2r] asit (2 — \ia.2 — j E ) [16]. In the following 
section we will apply the solutions (4.7a)-(4.7b) to specific inflationary models. We 
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will compare the accuracy of these solutions to that found by numerically integrating 
the full equations of motion for the perturbations. 

Before we proceed to the numerical results, we summarize the qualitative features 
of the solutions described here with typical time evolution of the slow-roll parameters 
found in §2. Using Eq. (2.10) and rjss ~ —0.1, the isocurvature perturbations grow 
steadily after the Hubble radius crossing and later source the curvature perturbations 
for 0(1/1^1) e-folds. In this case, the need for next-next-order contributions in Eqs. 
(4.7a)— (4.7b) comes from the fact that the expression —2r]^ s AN/3 can become a large 
fraction of r] ss if there are 0(l/\r)i^\) of e-folds between the Hubble radius crossing 
and the time when the coupling between the curvature and iscorvature perturbations 
reaches its maximal value. It also follows from Eqs. (4.5a)— (4.5c) that next-to-next 
order corrections should be negligible in this case. 

A different situation is encountered for the pattern described in Eq. (2.11). Here 
the isocurvature perturbations slowly decay due to nonzero 2e — r] aa until a spike in Re- 
produces a short increase in Vs which, as r} as is simultaneously driven to large values, 
can give a sizable sourcing of the curvature perturbations. This transient growth of 
the slow-roll parameters makes the use of the next order (or, generically, higher order) 
calculation necessary; whether this accuracy is sufficient depends on how large the slow- 
roll parameters actually become. As can be seen from Eqs. (4.5a)-(4.5c), the terms 
involving e come with large numerical coefficients, so when e becomes too close to 1, 
the splitting between different orders in slow-roll makes no sense and the perturbative 
calculation described here breaks down. 

5 Numerical examples 

Hybrid inflation [48] is unquestionably the most popular inflationary model with two 
scalar fields. In its original version, the inflaton slowly rolls down the potential until it 
reaches a critical value. Then, the waterfall field becomes tachyonic and quickly rolls 
down to the minimum of the potential, thereby terminating inflation. Such dynamics 
leads, however, to a blue spectral index, inconsistent with observations. This motivated 
a number of authors [33-35] to consider variants of hybrid inflation admitting a much 
milder transition between the inflaton-dominated and waterfall-dominated regimes of 
inflation 8 . In these models, the waterfall field remains tachyonic, but the magnitude of 
its mass parameter can be much smaller than the Hubble scale. In this case one has to 
analyze carefully the dynamics of the coupled system describing the behaviour of the 
perturbations of the two scalar fields. In §5.1 we investigate precisely these types of 
hybrid inflation models, whereas §5.2 will be devoted to models where the inflationary 

8 See [49] for a different approach where the authors investigate the importance of quantum diffusion 
on the transition. 
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Figure 1. The evolution of the slow-roll parameters e, rj aa , rj as and rj ss for the inflationary 
trajectory introduced in §5.1. 



trajectory has a turn in field space during which some slow-roll parameters may become 
close to one while maintaining the slow- roll condition e < 1. 

5.1 A case study of hybrid inflation 

As a concrete example, we shall study one particular inflationary trajectory which was 
also investigated in Ref. [35]. The potential of the model is 



V 2 J /i 2 ( fc v ' 2 



(5.1) 



where we adopt their notation and v = 0.10, ip c = 0.01 and /i = 1.00 x 10 3 . The overall 
normalization of the potential A does not play any role in this discussion and can be 
arbitrarily chosen to ensure the correct normalization of the curvature perturbations 
[3]. Our particular trajectory corresponds to the one shown in Fig. 4 of Ref. [35]: it 
starts at <p = 1.00 x 10~ 2 and Xo = 1-63 x 10 -9 , producing 62 e-folds of inflation. 
In the following, we shall mainly focus on the evolution of modes leaving the Hubble 
radius 8 e-folds after this initial time, so it is convenient to associate the beginning of 
the zeroth e-fold with that moment (negative number of e-folds will then correspond 
to earlier times). 

In Fig. 1 we show the time evolution of the slow-roll parameters e, r] aa , r] as and r) ss . 
It is clear that, for the initial 30 or so e-folds, the hierarchy (e, |r/ CTCT |, \f] as \) <C \rj ss \ < 1 
holds, realizing one of the possibilities discussed in §2.1. We note that the evolution 
of the slow-roll parameters r) aa , r\ as and rj ss qualitatively agrees with the approximate 
solution (2.10). 

We observe from Fig. 2 that the rf ss ^ corrections in Eq.(3.10c) are important for 
ensuring the accuracy of the result for a few e-folds. The results are normalized to the 
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Figure 2. Left panel: the power spectrum of the curvature perturbations, Vs, normalized 
to the single-field power spectrum of curvature perturbations V* calculated in four different 
ways as described in the text. Right panel: power spectra of the curvature and isocurvature 
perturbations (V-jz and Vs, respectively) for the case study of hybrid inflation described 
in §5.1. They were calculated numerically (red and yellow lines) or taken from analytical 
solutions (4.7a)-(4.7b) with parameters B and G expanded to leading (black dashed lines) 
or next-order (black solid lines) in the slow-roll parameters. The latter practically overlap 
with numerical results. 

value V* = H^/8ir 2 E+ of the single-field power spectrum of curvature perturbations. 
In the left panel, the evolution of Vs is calculated in four different ways. The yellow 
line corresponds to the numerical integration of the full equations of motion. The 
dashed black line is the prediction of Eq.(3.10c) neglecting rf ssi<1 while the solid black 
line corresponds to the full expression taking into account next-order effects. Finally, 
the thin grey solid line is the result of applying the transfer matrix formalism discussed 
in §3 with a leading-order slow-roll result for the power spectrum taken as the initial 
condition at N = 0. In these calculations, we have identified the number of e-folds N 
with ln(— fc*r), dropping the 0(e^) correction in Eq.(3.8). Since e is very small until the 
end of inflation (cf. Fig. 1), the error introduced by such an approximation is negligible. 
These results show that, for \r] ss ^\ = (9(10 _1 ), one must take into account corrections 
quadratic in i] ss in order to obtain a faithful estimate of the power spectrum of the 
isocurvature perturbations, valid for a couple of e-folds after Hubble radius crossing. 

Adhering to the particular trajectory in the model of hybrid inflation considered 
here, we have calculated the evolution of curvature and isocurvature perturbations 
during inflation for one particular mode leaving the Hubble radius at N = 0. The 
results are shown in Fig. 2 (right panel). Yellow and red lines correspond to the 
curvature and isocurvature perturbations, respectively, whose evolution was calculated 
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numerically from the full set of equations of motion (2.7a)-(2.7b). Black dashed (solid) 
lines show the predictions of our analytical solutions (4.7a)-(4.7b) with coefficients B 
and G expanded to next to leading order in the slow-roll parameters. 

We find good agreement between the full numerical solution and the approxi- 
mate solutions (4.7a)-(4.7b) obtained by including next-order slow-roll corrections. 
On the other hand, terminating the expansion at leading-order in slow-roll parame- 
ters gives inaccurate results. This behavior is a direct consequence of the fact that 
with r] ss ~ —0.1 the isocurvature perturbation is temporarily tachyonic, and grows 
for around twenty-five e-folds or so before the coupling between the curvature and the 
isocurvature perturbations becomes large, and the curvature perturbations are sourced 
by the amplified isocurvature modes. A correction of —2rf ss AN/2> in the exponent of 
(4.7b) can then easily be as large as around 20%, which is the magnitude of the effect 
we observe. 

In the 5N formalism, the curvature perturbations on large scales are given by 

n(t f ,x)=M(t f ,U,x)-N(tf,U) , (5-2) 

where Af(tf, £*, x) is the number of e-folds between an initial flat hypersurface at time t* 
and a uniform density hypersurface at time tf, and N(tf,t+) = J ' * H(t) dt. Following 
this definition, in two-field inflationary models one can write the power spectrum of 
the curvature perturbations as 

dN\ 2 _ fdN\ 2 



where and Q x are the Mukhanov- Sasaki variables associated with the perturbations 
of the scalar fields ip and \i respectively. Also, Vq and Vq x are the power spectra 
evaluated at and the partial derivatives are taken with respect to the values of the 
fields at t*. Eq. (5.3) does not rely on the slow-roll approximation. Usually, one further 
simplifies this result, associating t* with a time soon after Hubble radius crossing and 
taking Vq v = Vq x = H^/4ir 2 . However, for the inflationary trajectory considered 
here the latter assumption is not accurate enough: at Hubble radius crossing the 
power spectrum of the isocurvature perturbations is enhanced by nearly 20%, and this 
effect must be included in the calculation of the final power spectrum of the curvature 
perturbations. Upon identifying Q x with 5s, one notes that this correction is easily 
recognized as the 0(rf ss ) term in Eq.(3.10c). Importantly, this term should not be 
interpreted as a next-order correction to the 5N result, since the latter is already 
fully non-linear in the slow-roll expansion; instead it corresponds to the next-order 
correction to the initial conditions in the 5N formalism. 

To close this section we comment on how some of the subtleties discussed above 
impact existing results found in the literature. First, we note that Q a and 5s are 
two independent quantum fields. In order to obtain their power spectra, one should 
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solve the equations of motion for the wave functions twice, assuming either 5s = 
or Q a = as initial conditions; the two results should be then added in quadratures. 
Ref. [34] solves the equations of motion for the wave functions only once, which leads 
to a spurious dip in the plot of \Q a \ shown in logarithmic scale. However, thanks 
to a fortunate coincidence, the final result is numerically correct; the final curvature 
perturbations are dominated by the contribution of the isocurvature perturbations, so 
the initial curvature perturbation can be safely neglected. In Ref. [35] the sourcing 
of the curvature perturbations by the isocurvature perturbations was neglected, which 
not only leads to a final amplitude of the curvature perturbations that is orders of 
magnitude smaller than the full result, but also gives an incorrect estimate n s — 1 = 
— 6e* + 2r) aa * for the scalar spectral index. It follows from (4.7a) that [16] 

n s - 1 ~ -(6 - 4c 2 K + 2(1 - c 2 )7] aa , + 2c V* , (5-4) 

where c 2 = 1 — V^/V-ji is the isocurvature-sourced fraction of the curvature perturba- 
tions (here we assumed c* = 0). However, for the trajectory discussed here c 2 ~ 0.999, 
so n s rH 0.8, which is excluded by WMAP results [3] (irrespective of the constraints 
from the normalization of the perturbations [50]). More generally, with negligible £*, 
rj uai , and T] asin we find that c 2 = (n s — l)/2r] 89if , so with r] ss * ~ —0.1 and n s — 1 ~ —0.04 
we can only have c 2 ~ 0.2. Fig. 2 shows that, compared to the leading-order results, 
for 7] s s * ~ — 0.1 we have ~ 10% correction to c 2 , coming from a similar correction to 
Vs*- From this estimate we can conclude that for phenomenologically viable models 
with the pattern of the slow-roll parameters considered here, the next-order corrections 
at Hubble radius crossing result in corrections to the power spectrum of the curvature 
perturbations at the level of a few percent. This is comparable to the expected accuracy 
of the Planck measurements. 

5.2 A case study of double quadratic inflation 

Double quadratic inflation models [44], given by the potential 

V(^,X) = ^ml^> 2 + ^m 2 x x 2 , (5.5) 

are a well-studied class of multi-field inflation models. Here we study a trajectory 
closely resembling the example investigated in Ref. [17]: we take m x /m Lp = 7 and 
start the evolution of the homogeneous fields at ipo = 8 and xo = 8 (the overall 
normalization of the potential only affects the normalization of the perturbations). 
Inflation is initially driven by \- Then, around N = 10 the field space trajectory 
rapidly changes direction and there is a slight deviation from slow-roll, as measured 
by the e parameter 9 . 

9 In the models considered here and in §5.3, the number of e-folds is smaller than required by 
observations. However, these models are only intended as illustrations of the general effects we are 
studying. 
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Figure 3. Case study of double quadratic inflation described in Sec. 5.2. The left panel 
shows the evolution of the slow- roll parameters e, rj aa , tj^s and r/ ss . The right panel shows 
the power spectra of the curvature and isocurvature perturbations (V-r and Vs , respectively) 
calculated numerically (red and yellow lines) or taken from the analytical solutions (4.7a)- 
(4.7b) with parameters B and G expanded to leading (black dashed lines) or next-order 
(black solid lines) in the slow-roll parameters. 

The evolution of the slow-roll parameters e, r] aa ., r\ as and r\ sa is shown in Fig. 3, 
where we also display the evolution of the perturbations. Since in this example 
i"j JK = 0, the coupling B between the curvature and isocurvature perturbations on 
super-Hubble scales, given by Eq.(4.4b), is proportional to r] as , which explains a small 
decrease in Vn after reaching the maximum at N ~ 11. The slow- roll expansion be- 
comes unreliable when r] ss becomes much larger than one. Our next-order slow-roll 
expansion is, however, quite reliable before that happens. The discrepancy between 
the results calculated up to leading and next-order in slow-roll comes from the con- 
tributions to the parameter B, with the next-order contribution only twice smaller 
than the first order one. At later times the isocurvature perturbations become very 
massive and decay rapidly, no longer affecting the evolution of the curvature modes. 
The most striking result of the breakdown of the slow-roll expansion is the oscillatory 
feature in the power spectrum of the isocurvature perturbations Vs, resulting from a 
sudden change of the mass parameter 3rj ss H 2 of this perturbation. This also explains 
an interesting feature in Fig. 3: after decaying for a short period (when N ~ 11) the 
quantity Vn increases slightly. This effect goes beyond our analytic approximation, 
but has almost no impact on the final value of the power spectrum. 
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Figure 4. Case study double quartic inflation described in Sec. 5.3. Left panel shows 
the evolution of the slow- roll parameters e, r) aa , r} as and r] ss . Right panel shows the power 
spectra of the curvature and isocurvature perturbations (Viz and Vs, respectively) calcu- 
lated numerically (red and yellow lines) or taken from analytical solutions (4.7a)-(4.7b) with 
parameters B and G expanded to the leading order (black dashed lines), next order (black 
solid lines) or next-to-next order (black dash-dotted lines) in the slow-roll parameters. 

5.3 A case study of double quartic inflation 

This example originates from the model of matrix inflation [51]. With certain simpli- 
fying assumptions the potential of this model can be written as 

V( V , X )= 1 1 \^ + \\ X X 4 - (5.6) 

We investigated the second of the exemplary inflationary trajectories studied in [51], 
taking \ ip /\ x = 410 with ipo = 11.2 and xo — 9.1. Initially, inflation is driven by tp, 
then, around N = 8, the trajectory in the fields space rapidly changes direction and 
there is a deviation from slow-roll, as measured by the e parameter. 

The evolution of the slow-roll parameters e, r] aa , r] as and r] ss is shown in Fig. 4, 
where we also display the evolution of the perturbations. It is a clear example of the 
limitation of the method presented in §4: the higher order terms in Eqs. (4.5a)-(4.5c) 
contribute terms with coefficients much larger than one and containing powers of e 
in Eqs. (4.7a)-(4.7b). Hence whenever e becomes too close to 1 during inflation, the 
expansion in the power series of the slow-roll parameters starts to break down. The 
example shown here is on the verge of applicability of the perturbative expansion in 
the powers of the slow-roll parameters, and it can be seen from Fig. 4 that the next-to- 
next order provides a fairly improved estimate of the curvature and the isocurvature 
perturbations. 
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6 Conclusions and outlook 



In the present era of precision cosmology, the Planck satellite should deliver high- 
precision data in the near future. Such observational standards need to be followed 
by equally good theoretical precision so that one can make the most efficient use of 
the Planck data. Whilst the inflationary power spectrum at lowest-order in slow-roll 
can provide an accurate description of cosmological evolution in many scenarios, there 
exist models where one requires slow-roll effects beyond lowest-order for an accurate 
description of the evolution of the perturbations. 

In this paper we have studied the evolution of the curvature and isocurvature 
perturbations in two-field inflation models to next-next-order in the appropriate slow- 
roll parameters. Using popular examples of two-field inflationary models, a variant of 
hybrid inflation and double quadratic/quartic inflation, we have found that calculating 
the two-point correlations of the isocurvature modes at next-order in r) ss provides a very 
precise estimate of the initial amplitude of the perturbations around the Hubble radius. 
The isocurvature perturbations can source the curvature perturbations during a turn 
in the inflationary trajectory in field space: all slow-roll parameters can (temporarily) 
increase at this instance, which justifies keeping next-order effects in the analysis. To 
follow the dynamical evolution, we have applied the transfer matrix method and have 
found a remarkable agreement with the explicit numerical integration of the equations 
of motion. 

In conclusion, our results show that keeping slow-roll corrections to higher or- 
der is necessary to obtain an accurate description of perturbations on super-Hubble 
scales. We reached this conclusion using knowledge of a single inflationary trajectory, 
which had the advantage of making the precision of our semi-analytical calculations 
comparable to those of the full numerical methods. In the examples we considered we 
found that this enables this class of models to be predictive at an accuracy required 
for comparison with data expected from Planck. 
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A Equations of motion for the perturbations 

Some of the results reviewed in §2 and §3 require rather lengthy expressions which 
might obscure the main line of thought. For completeness of the discussion, we have 
collected them here. 

The coefficients Cjj in the equations of motion (2.7a) and (2.7b) are: 

C„ = V^-(^) + 2a^ + 3& 2 -^-, (A.la) 



a 



H 2H 2 



C as = 6H^ + 2^4^ + 2V as + , (A.lb) 
(J <J ri 



C sa = -6H^ - 2 V -^f + , (A.lc) 
cr cr H 



v- 2 



c ss = v ss - > ( A - ld ) 



a 

which are exact to all orders in slow-roll. We can expand these coefficients uniformly 
to next-order in the slow-roll parameters: 

Caa^ 3H 2 (^ aa - 2e - l -r,l s + ^e Vaa - ^ , (A.2a) 

C as ~ \2r] as + ^VasVaaj , (A.2b) 

C sa ~ 3H 2 (-^r) as r) a(T + ^ £ Vas) > (A.2c) 
C ss ~ 3H 2 Us - \ril)j • (A.2d) 

We can proceed similarly with respect to the coefficient V s /a, multiplying the first 
derivatives of the perturbations as Hr} as {l + e/3). 

One can now remove the first derivatives from the system of equations (2.7a) and 
(2.7b) by performing a change of basis. We introduce the rotated field perturbations 

Wl )=R-i( U °) (A.3) 
w 2 J \u s J 

where the rotation matrix is defined in terms of an angle d 

R-( C0S ! ~**?) , (A.4) 

\ SUlt/ COS 17 / 
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with 

■&=-£. (A.5) 
a 

The equations of motion for the perturbations w\ and u> 2 around the time of horizon 
crossing are, to next-order in the slow-roll parameters 



+ 



k 2 - -) 1 + -Ir^qr 



a / t 2 



Wi 

w 2 



0, (A.6) 



q _ I 3?7 (TCT - 6e + 10er] aa 3r] as + 8er] as . 



where 

/ X'jn f\ cr _1_ 1 (icrrt <in _1_ era , 

The assumed pattern of the slow-roll parameters justifies treating them as constant in 
the course of a few e-folds, around the time when the mode of interest with a comoving 
wave number crosses the Hubble radius, that is, k* = a*H*. We can therefore replace 
the slow-roll quantities by their values at Hubble radius crossing. Moreover, •d can also 
be replaced by its value d± at Hubble radius crossing, and we have the freedom of 
choosing the constant term in the solution of (A.5) so that the matrix R _1 QR is 
diagonal. We can then write the equations of motion for the perturbations W\ and w 2 
as 

/ II \ T -t / t-* , \ r\ \ ~l / 

W 1 



// 

2 



W 



+ 



-2 



2 + A 2 * 



W'2 



, (A.* 



where we used a" /a = (2 — e)(aH) 2 . Also, when writing Eq.(3.2), we used the exact 
relation fact V s /& = Hrj as /(1 — e/3). The linear combinations 

A A * = 3e* + 20el + 8e^ ffM - A^* (A.9) 

where A e {1, 2}, can be expressed in terms of the angle as 

(Ai* - A 2 *) sin(2$*) = 677^ + 166*77^ , (A. 10a) 

(Ai* - A 2 *) cos(2#*) = 37/^ - 3t7 ss * - 6e* + lOe+r)^ + 6£*r] ss ± . (A. 10b) 

The relations 

Ai* + A 2 * = 3ri aa + + 3r/ ss * - 6e± + 10e*r/ CTCT7t + 6£*T7 SS * , (A. 11a) 

AL + AL ^ HvL + 2r?L + Via* ~ 4e*r)a<r* + 4e^) , (A.llb) 

are also useful for obtaining the two-point correlations of the curvature and isocurva- 
ture modes. With the assumed hierarchy e, \r) aa \, \r) as \ <C \r) ss \ < 1 the results above 
can be further simplified: 

Ai = 3r] aa ^ - 6e* + lOe^^ - 9^ , (A. 12a) 

A 2 = 3r/ SSVc + 6e*r/ ss * . (A. 12b) 
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Note that all the previous equations have been expanded uniformly to next-order in 
the slow-roll parameters. However, for the case studied in §5.1, only terms quadratic in 
the largest slow-roll parameter, rj ss , are relevant for maintaining the required accuracy 
in the evolution of perturbations. The only quadratic contributions to the equations 
of motion (3.2) arise from quadratic terms in A^. The case studied in §5.2 has all 
the slow-roll parameters much smaller than 1 at Hubble radius crossing, therefore the 
standard results apply. 
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